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Abstract 

A large luminosity-linewidth template sample is now available, improved 
absorption corrections have been derived, and there are a statistically sig- 
nificant number of galaxies with well determined distances to supply the 
zero point. A revised estimate of the Hubble Constant is Ho = 77 ± 4 km 
s _1 Mpc^ 1 where the error is the 95% probable statistical error. Systematic 
uncertainties are potentially twice as large. 

1. Introduction 

More massive galaxies have more stars than less massive galaxies and more 
massive galaxies rotate faster. This simple reasoning leads to the expecta- 
tion of a correlation between the light of the stars and a measure of the 
rotation rate (Tully & Fisher 1977) though the small scatter of the cor- 
relation is not trivially explained. The measurement of the rotation rate, 
say from the width of a global neutral hydrogen profile, is independent 
of distance. Hence, if the relationship is calibrated in terms of the intrin- 
sic luminosity dependence with linewidth then the modulus between the 
apparent and absolute magnitude at a given linewidth gives a distance. 

Nowadays there are more accurate methods for measuring the distance 
to an individual galaxy. However, the great virtue of the luminosity-linewidth 
method is that it can be used to get distances for thousands of galaxies 
all over the sky. Roughly 40% of galaxies with Mb < — 16 m are poten- 
tial targets. Current techniques allow candidates to be accessed out to 
V ~ 10,000 km s _1 . Consequently, a determination can be made of the 
Hubble Constant based on distance measurements dispersed on the sky 
and in a regime with modest consequences due to peculiar velocities. 
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It is a good moment to review the results obtained with this method. 
Lately there has been a big input of high quality data. Coupled with this 
happy situation is the dramatic improvement in the absolute calibration as 
a result of the determination of distances by the cepheid pulsation method 
with Hubble Space Telescope. These two improvements address what were 
the greatest deficiencies in the path to the Hubble Constant based on the 
luminosity-linewidth method. It will be shown during the ensuing discus- 
sion that today there are sufficiently large numbers of both calibrators and 
targets that the statistical accuracy of the method is very good. Uncertain- 
ties are now dominated by potential systematic effects. 

Here is an outline of the structure of the article. First there will be 
a discussion of the raw data: luminosities, axial ratios, and linewidths. 
The information comes from many sources. The observed parameters re- 
quire adjustments for modifying effects like extinction and projection. Then 
the methodology of the construction and application of the luminosity- 
linewidth correlation will be described. The potential problems of Malmquist 
bias must be confronted. It will be described how a template luminosity- 
linewidth correlation is constructed and then transformed to an absolute 
magnitude scale. Then the template relation is imposed upon data in clus- 
ters that for the most part are beyond the Local Supercluster, at distances 
where peculiar velocities should be only a modest fraction of expansion ve- 
locities. That exercise results in distance estimates for the clusters and a 
determination of the Hubble Constant. 

2. Data 

Three parameters must be measured: an apparent magnitude, a character- 
ization of the rotation rate, and an estimator of the inclination needed to 
compensate for projection effects. Each of these components will be con- 
sidered in turn. Then, there will be a discussion of the adjustments to be 
made to get to the parameters that are used in the correlations. 

2.1. LUMINOSITIES 

Area photometry with optical and near-infrared imagers has come of age. 
Large format detectors on modest sized telescopes provide fields of view 
that encompass the entire target galaxies. The author has had an on-going 
program of both optical and infrared photometry (Pierce & Tully 1988, 
1992; Tully et al. 1996, 1998). For the purposes of the present discussion, the 
other important sources of luminosities are Mathewson, Ford, & Buchhorn 
(1992), Han (1992) and Giovanelli et al. (19976). The latter three sources 
provide / band magnitudes for galaxies in clusters at intermediate to large 
distances. The collaboration involving Pierce and Tully is producing B,R,I 
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magnitudes for many nearby galaxies, including the calibrators, and for 
galaxies in a couple of distant clusters that overlap with the other sources. 
In addition, data in the K' band is provided in Tully et al. (1996, 1998) for 
two of the clusters. 

At the moment, there is a lot more data available at / than at other 
bands so most of the analysis presented in this paper will be based on 
this material. There is interest in the other bands, though, because of the 
insidious effects of obscuration. There will be some comfort that there is 
proper compensation for these effects if relative distances are the same at 
different passbands. The K' material is of particular interest in this regard 
since obscuration should be very small at 2 microns. 

The issue of adjustments to magnitudes because of obscuration and 
spectral shifting will be discussed in a later section. The concern at this 
point is the homogeneity of the raw magnitudes from various sources. Dif- 
ferent authors measure magnitudes to slightly different isophotal levels then 
usually extrapolate to total magnitudes: Han extrapolates from 23. 5 m , Gio- 
vanelli et al. extrapolate from ~ 24 m , Mathewson et al. extrapolate from 
25.0 m , Tully et al. (1996) extrapolate from 25. 5 m , and Pierce & Tully give 
a total magnitude to sky at ~ 26 m . The added light at the faintest levels is 
small for the high surface brightness galaxies that are relevant for the deter- 
mination of Ho. Typical extrapolations from 25. 5 m to infinity add ~ 0.02 m 
and always less than 0.1 m for galaxies in the appropriate magnitude range 
(Tully et al. 1996). Magnitudes measurements are more vulnerable to the 
detailed fitting of the sky level. Variations at the level of ~ 0.05 m can arise 
with systematic differences in sky fitting procedures. 

Inter-comparisons between sources indicate that everybody is working 
on the same system and that systematics are almost negligible. Some offsets 
have been reported, for example Giovanelli et al (19976) adjust Mathewson 
et al. data (1992) to match their own. However, globally the data sets are 
consistent with each other at a level of 2% in effect on distances. Object by 
object, rms differences between any pair of observers is at or below ±0.1 m . 
In the present analysis all sources are given equal weight and luminosi- 
ties are averaged if there are multiple observations. Overlap measurements 
reveal spurious results in a few percent of cases. If a difference between 
sources is big enough it is usually evident which measurement is wrong. 

2.2. INCLINATIONS 

Projection corrections are required to recover true disk rotation rates and 
to compensate for differential obscuration. Uncertainties in inclination es- 
pecially affect de-projected velocities as one approaches face-on. With rare 
exception, inclinations are derived from a characteristic axial ratio of the 
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main or outer body of a galaxy. From experience, it is found that such 
inclination measurements are reproduceable at the level of ±3° rms. How- 
ever, errors are non-gaussian. From the radial variations in axial ratios and 
from such independent considerations as inclination estimates from velocity 
fields it is suspected that errors of ~ 10° are not uncommon. The 1/sini de- 
projection correction becomes grossly uncertain toward face-on. A sample 
cut-off at i = 45° is invoked to avoid large errors. 

The derivation of an inclination from an axial ratio requires an assump- 
tion about the intrinsic thickness of the system. The standard formulation 

for the derivation is cosi = \J (g 2 — g 2 )/(l ~~ Qo) where q = b/a is the ob- 
served ratio of the minor to major axes and qo is the intrinsic axial ratio. 
The thinnest systems are spirals of type Sc. Earlier types have bulges and 
later types are puffed up. For simplicity, go = 0.20 is often used. A more 
elaborate specification of qo that depends on type could be justified. Gio- 
vanclli et al. (19976) provide an extreme example with their choice qo = 0.13 
for type Sc. A smaller qo value results in derived inclinations that are more 
face-on. Fortunately the choice of qo has negligible effect on the measure- 
ment of distances as long as one is consistent. The difference go = 0.13 or 
0.20 gives a difference in inclination, for an observed q = 0.20, of 81° or 
90° respectively. However the 1 / sini difference on the corrected linewidth is 
only 1.2%. As one progresses toward larger q the difference in assigned in- 
clination is reduced but the 1/sim amplifier is growing. The product of the 
two effects is a roughly constant shift of 1.2% in the corrected linewidth at 
all inclinations i > 45° . If both calibrators and distance targets are handled 
in the same manner there will be no effect on measured distances. 

Difficulties with projection enter luminosities in the opposite regime, as 
galaxies are presented toward edge-on. It has become popular to formulate 
extinction corrections directly in terms of the observed q value which avoids 
a dependence on the parameter go. 

Inter-comparisons between the sources of photometry used in this study 
fail to reveal any systematic differences in g measurements between authors, 
though big individual differences are not uncommon. Big differences raise 
flags that prompt special attention. 

2.3. LINEWIDTHS 

It has become popular to measure rotation parameters via both optical 
and radio techniques. The original radio methods are simpler but are con- 
strained by detector sensitivity to modest redshifts. The methods that in- 
volve optical spectra require more work but can be used to larger distances. 
With care, the two techniques can be reconciled in a common characteriza- 
tion of the projected rotation. However that synthesis will not be attempted 
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here. There are plentiful and sufficiently distant observations of profiles in 
the 21 cm neutral hydrogen line for the purpose of determining Ho- The 
complexity of intermingling radio and optical data can be avoided. 

Even restricted to the HI data, there is a bit of a mess. For once in 
astronomy, resolution is not an unmitigated advantage. In this case it is 
desired that the beam project onto an area larger than the galaxy in order 
to enclose most of the emission. As a result, the data on nearby, large 
galaxies is handed down from observations on old telescopes from the days 
of paper strips. The parameterizations are still quite 'personalized'. Worse 
than with magnitudes, one has to be careful that one is using a consistent set 
of linewidth information from the near field to far. In this study, HI profile 
linewidths defined at the level of 20% of the peak flux are used (called 
W20). These linewidths are only adequately measured if the line signal-to- 
noise is greater than 7. A decent signal typically provides a measurement 
with an accuracy of better than 10 km s _1 . The 20% linewidths are then 
transformed into the parameter Wr defined by Tully & Fouque (1985). This 
parameter approximates twice the maximum rotation velocity of a galaxy. 

The other linewidth characterization in common use is the width at 50% 
of peak flux in each horn of the profile (W50: Haynes et al. 1997) which is 
then adjusted to account for instrumental and thermal broadening (Gio- 
vanelli et al. 19976). The advantages and disadvantages of the alternative 
systems are technical and not important. The key concern is that the in- 
formation available over both north and south hemispheres and for both 
nearby large galaxies and those distant and small be brought to a common 
system. The current analysis draws on a large database of W20 measure- 
ments and a supplement of W50 values in clusters well beyond the Local 
Supercluster. The important new contributions in the W50 system are only 
partially within the public domain so it has not yet been possible to make a 
detailed comparison based on the ~ 10 3 galaxies that must be mutually ob- 
served in the two systems. The present study relys on an inter-comparison 
of only 66 galaxies in 3 clusters. It is found that < W20 — W50 >= 25 km s _1 
with 12 km s _1 rms scatter and a standard deviation of 2 km s -1 . There is 
no apparent trend with W or between the three clusters. While acceptable 
for a preliminary foray, this inter-comparison of the W20 and W50 systems 
must, and can easily, be improved upon. 

2.4. EXTINCTION CORRECTIONS 

Along with the Hubble Space Telescope contribution to the zero-point cal- 
ibration and the abundance of new material, the third significant improve- 
ment of late has been in the compensation for extinction. Giovanelli et al. 
(1995) made a convincing case for a strong luminosity dependence in the 
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Figure 1. Dependence of the extinction amplitude parameter 7a on absolute magnitude. 
Data are presented for B,R,I bands in the three separate panels. The filled squares 
correspond to values of 7a derived from deviations from mean color relations as a function 
of b/a. The filled triangles correspond to equivalent information derived from deviations 
from luminosity-line profile width correlations as a function of b/a. The small circles in 
the I panel are data taken from Giovanelli et al. (1995). The dashed straight line in the 
I panel is a least squares fit to the Giovanelli et al. data (errors in 7). The solid straight 
lines in this and the other panels are least squares fits to the data by Tully et al. (1998). 
The dotted straight line in the / panel gives equal weight to the two surces of data. 



obscuration properties of galaxies and Tully et al. (1998) have further quan- 
tified the effect. The latter work has profited from the leverage provided 
by information in passbands from B to K' . A giant galaxy can be dimmed 
by 75% at B if it is viewed edge-on rather than face-on, although a dwarf 
galaxy of the luminosity of the Small Magellanic Cloud statistically has 
not enough extinction at B to measure. At K' the most luminous galaxy 
is dimmed by a maximum of 20%. 

The extinction can be described by the expression A\ = 7\log(a/6) 
where a/b is the major to minor axis ratio and A is the passband. The 
correction is to face-on orientation but does not account for the residual 
absorption in a face-on system. Figure 1 provides a plot from Tully et al. 
(1998) that shows the dependence on luminosity of 7a for A = B,R,I. 
Given this strong luminosity dependence, there is a problem because abso- 
lute magnitudes are not known a priori. Absolute magnitudes are to be an 
output of the distance estimation process so they cannot also be an input. 
Both Giovanelli et al. (19976) and Tully et al. (1998) recast the corrections 
for magnitudes so the dependency is on the distance independent linewidth 
parameter. This conversion is provided through the luminosity-linewidth 
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calibrators. The formulations presented by Tully et al. (1998) are: 

7 B = 1.57 + 2.75(logW£ - 2.5) (1) 

1R = 1.15 + 1.88(logWj - 2.5) (2) 

7/ = 0.92 + 1.63(logW& - 2.5) (3) 

lK , = 0.22 + 0.40(logWjj - 2.5) (4) 



There is a fortunate interplay that minimizes the effect of uncertain inclina- 
tion on A\. If the inclination is taken too face-on because of an spuriously 
large b/a then W l R is overestimated, which drives up 7a, but is offset by a 
low log (a/6) in the product that gives A\. 

The other corrections to be made are modest and non-controversial. 
Absorption at I due to obscuration in our own Galaxy is taken to be 41% 
of the B band value given by Burstein & Heiles (1984). There is a small 
'k-correction' of 1.27.Z. 

3. Methodology 

Over the years many people have used luminosity-linewidth relations to 
measure distances and there has been controversy. An extreme view has 
been presented by Sandage (19946). According to him, there can be large 
biases that distort distance measurements and limit the usefulness of the 
procedure. In this section there will be a description of a way of conducting 
the analysis that results in unbiased distance estimates and, hopefully, ac- 
curate results. The reader interested in making a comparison will find that 
the method to be described is not the method used by Sandage. 

3.1. BIASES 

Malmquist (1920) discussed a bias that might create a problem with mea- 
surements of distances to objects selected by apparent magnitude. Teeriko- 
rpi (1984) and Willick (1994) have discussed the problem in the present 
context. Schechter (1980) and Tully (1988a) have described a procedure 
that is expected to nullify the bias. That procedure will be summarized. 

An example of when the bias arises is provided by considering the de- 
scription of the luminosity-linewidth correlation given by the regression 
with errors taken in magnitudes - sometimes called the 'direct' relation. 
Use the 'direct' relation to determine distances to objects in the field. By 
the construction of the regression, the brightest galaxies will tend to lie 
above the correlation line. Suppose one considers a group. The brightest 
galaxies, drawn from above the mean correlation but assigned the absolute 
magnitude of the mean correlation, will be given a closer distance than 
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is correct. As fainter galaxies in the group are sampled, they progressively 
sample the true distribution around the mean correlation, so that the mean 
distances of the fainter galaxies are larger. Kraan-Korteweg, Cameron, &; 
Tammann (1988) have shown that the measured mean distance of a group 
increases as fainter objects are included. For the same reason, as one probes 
in the field to larger redshifts one samples progressively only the brightest 
galaxies, those that tend to be drawn from above the mean correlation. 
Hence one progressively assigns erroneously low distances. Low distances 
give a high Ho- 

In an analysis made this way it is imperative that a correction be made 
for the bias. However, to make the correction it is necessary to have detailed 
information on the form of the luminosity-linewidth correlation and the 
nature of the scatter. With adequate information, it is possible to correct 
statistically for the bias, though the trend of deviations with magnitude 
would persist in the individual measurements. It is submitted that Sandage 
(19946) provides an example where the characteristics of the correlation and 
scatter are not understood and the corrections are erroneous. 

Variations on the procedures that require bias corrections are pervasive 
(eg, Willick et al. 1997). For example, a maximum likelihood description of 
the relationship (Giovanelli et al. 19976) still retains the bias and requires 
corrections. The corrections might be done properly. However, these proce- 
dures require (1) that the calibrators and targets have the same statistical 
properties, and (2) detailed specification of the sources of scatter and of 
properties of the luminosity function from which the sample is drawn. As 
an alternative, the method to be described nulls the bias rather than cor- 
rects for it. Consequently, there is no requirement to specify the sources of 
scatter or the properties of the sample. One is relying only on the assump- 
tion that calibrators and targets have the same properties. 

The magic description that nulls the bias is given by the regression with 
errors in linewidth (Schechter 1980; Tully, 1988a, 6) - the 'inverse' relation. 
Two qualitative comments might crystallize the merits of the procedure. 
The first point to appreciate is that the amplitude of the bias depends on 
the assumed slope of the correlation. The flatter the dependence of magni- 
tude with linewidth the greater the bias. Conversely, if the slope is taken 
steep enough the sign of the bias can be reversed. Hence it can be under- 
stood that there is a slope that nulls the bias. That slope is given by the 
regression on linewidth if the sample is only limited in magnitude. The 
second key point is made by a consideration of the regressions on the sepa- 
rate axes of a luminosity-linewidth plot. Suppose one considers successively 
brighter magnitude cuts on an intrinsic distribution. As one progressively 
limits the magnitude range, the correlation coefficient of the fit will degrade. 
Presented graphically, the correlations on the two axes will progressively 
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diverge as the fitting range is reduced. Here is the critical point. As the 
truncation is progressively advanced in magnitude the slope of the correla- 
tion with errors in linewidths is always the same but the slope with errors 
in magnitudes is progressively splayed to shallower values. 

Since the amplitude of the bias depends on the slope of the correlation, 
it should be seen that an analysis based on the direct relation is on slip- 
pery ground because the value of the slope depends on the magnitude limit 
of the sample. One needs a lot of information for an internally consistent 
application. The maximum likelihood approach raises the same qualitative 
concerns although, because it involves a slope intermediate between the di- 
rect and inverse correlations, the quantitative problem is also intermediate. 

It has been pointed out by Willick (1994) that a bias can enter the 
inverse correlation in practical applications. The bias can be introduced 
because the cutoff may not be strictly in magnitude. For example, the 
sample might be chosen at B band but applied at a more redward band 
such as /. A correlation between color and linewidth generates a slope to 
the magnitude cutoff at a band other than B. Or suppose the sample is 
selected by apparent diameter. A correlation between surface brightness 
and linewidth can again give a slope to the magnitude cutoff. A slope 
in the magnitude cutoff is equivalent to the introduction of a linewidth 
stricture. Any restriction in linewidths brings the problem of bias over to 
the orthogonal axis. Two things can be said of this problem. First it is a 
small effect, down by a factor of five in amplitude in Willick's analysis. 
Second the problem is partially avoided by building the calibration out of 
only galaxies that satisfy a completion limit at the band to be considered; 
ie, a stricter limit is taken than the one that provided the initial sample. 

Most important: to achieve the correlation that nulls the bias one wants 
a complete magnitude limited calibration sample. In the population of the 
luminosity-linewidth diagram with the calibration sample there should not 
be any discrimination against candidates in any particular part of the dia- 
gram above the magnitude limit. Selection based on inclination is inevitable 
but that restriction should be distributed across the diagram. Other poten- 
tial restrictions must be considered in a similar light. 

The good news is that, with due care to the calibration, then the method 
can be applied to give unbiased distances to individual galaxies in the field 
as long as the inclusion of those galaxies is not restricted in linewidth. In 
other words, there will not be a correlation between luminosity and distance 
within a group as found be Kraan-Korteweg et al. (1988) nor a correlation 
between Ho and redshift as found by Sandage (1994a). The method will 
break down if the target galaxy is a dwarf intrinsically fainter than the 
limit of the calibration. The latter issue is only a concern in our immediate 
neighborhood, not for the Hq problem. 
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3.2. THE TEMPLATE RELATION AT I BAND 

The creation of the template relation is a critical step. In the section on 
biases it has been described how important it is to have a sample that only 
suffers magnitude constraints. Often the calibration relationship is formed 
out of the ensemble of a field sample (Willick et al. 1996) but the constraints 
on such samples are usually ambiguous. Also, the calibration relationship 
is inevitably broadened and distorted by non-Hubble expansion motions. 

Cluster samples have evident advantages. It is possible to be complete 
to a magnitude limit and it can be assumed that the galaxies are all at the 
same relative distance. The biggest concern with cluster samples is whether 
there are intrinsic differences between galaxies in a cluster environment and 
those that are more isolated. An operational disadvantage of cluster samples 
is that an individual cluster does not provide enough systems to provide 
good statistics. These two disadvantages can be addressed simultaneously 
by building a template relation out of several cluster samples. The 'clusters' 
can have a sufficient range in their properties that one can begin to evaluate 
the issue of environmental dependence. The combination of several cluster 
samples takes care of the problem of poor statistics. 

This study uses samples drawn from five clusters with reasonable com- 
pletion characteristics. There is best control with the nearby Ursa Major 
and Fornax clusters. The completeness limits in Ursa Major are discussed 
by Tully et al. (1996) and in Fornax by Bureau, Mould, and Staveley-Smith 
(1996). After corrections for obscuration, and translation to I magnitudes, 
the completion limit for both clusters is / = 13. 4 m . There are 38 galaxies 
in Ursa Major with type Sa or later and i > 45° above this limit. There 
are 16 galaxies in Fornax satisfying these constraints. It was appreciated in 
advance that Ursa Major and Fornax are at similar distances. Hence the 
apparent magnitude limits conform to about the same absolute magnitude 
limits. Fornax is indicated by these data to be 0.10 m closer. 

Already a diverse environmental range has been explored between the 
Ursa Major and Fornax cases. Tully at al. (1996) have labored the point 
that the Ursa Major Cluster environment is more similar to that of low 
density spiral groups than to what is generally considered a cluster. The 
structure must be dynamically young. By contrast, Fornax has a dense 
core of early type systems, evidence of a dynamically evolved structure. 
Granted, the spirals in the Fornax sample are more widely distributed than 
the central core and may represent recent arrivals. 

The next component to be added to the template is drawn from the fila- 
ment that passes through what has been called the Pisces Cluster. Aaronson 
et al. (1986) and Han h Mould (1992) have included the region in their dis- 
tance studies but Sakai, Giovanelli, <fe Wegner (1994) have shown that one 
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is dealing with an extended structure with separate sub-condensations. It is 
unlikely that the region as a whole is collapsed. Indeed, what will be consid- 
ered here is a length of ~ 20° along the Pisces filament, which corresponds 
to an end-to-end distance of ~ 20 Mpc. The mean redshift is constant to 
~ 4% along the filament though individual redshifts scatter over a range 
of ±20% relative to the mean. It can be asked if the full length of the fila- 
ment is at a common distance or if variations in distance can be identified. 
A luminosity-linewidth correlation is constructed for the ensemble, then 
inter-compared by parts to determine if components deviate significantly 
from the mean. There is not the slightest hint of deviations from the mean. 
Six sub-components along the 20° filament have consistent distances within 
a few percent. To within measurement errors, the filament is tangent to the 
plane of the sky in both real space and velocity space. 

Given this circumstance, all the galaxies with 3700 < V cm b < 5800 km s~ 
along the 20° segment of the Pisces filament 00 h 44 m < a < 02 fe 13 m will 
be taken to be at the same distance. Failures of this assumption will act to 
increase the scatter of the luminosity-linewidth relationship but the scat- 
ter is found to be only 0.31 m , as small as for any sub-component of the 
template. This scatter is with 46 galaxies, after rejection of one object that 
deviates by ~ 4a. There is reasonable completion brighter than / = 13. 8 m 
which is taken as the magnitude limit for the present sample. The Pisces 
filament data is added to the Ursa Major /Fornax template by (1) calculat- 
ing the offset from the slope of the 2 cluster template, (2) redetermining 
a new slope now with 3 clusters, (3) iterating the distance offset with the 
new slope, and (4) iterating the new 3 cluster template slope. The distance 
shift at step 3 is of order 1% and the slope shift at step 4 is ~ 1%. 

The final step in the development of the template is the addition of the 
Coma and Abell 1367 clusters. These clusters are at the same distance to 
within a few percent so they are treated together until the final iteration, at 
which point they are considered separately against the mean relation. Only 
galaxies within 4.3° of the cluster centers are accepted and the velocity 
constraints described by Giovanelli et al. (19976) are accepted. As with 
Pisces, there is substantial but not full completion to / = 13. 8 m . Iterations 
like those described with the Pisces filament converge to provide the five 
cluster template. There could have been a problem if there is curvature in 
the template, as might be indicated if, say, the slope flattened for samples 
with more luminous cutoffs (more distant clusters). However there is no 
suggestion of such a flattening. The Coma sample provides 28 galaxies. 
A1367 adds 23, after one 5a rejection. 

In total, there are 151 galaxies in the 5 cluster template. There are three 
distinct absolute magnitude cutoffs (UMa/Fornax; Pisces; Coma/A1367) 
but, to the degree that the slopes are indeed constant between components, 
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Figure 2. Template luminosity-linewidth relation at I constructed with 151 galaxies in 
five clusters, translated in zero-point to get a best fit with 16 calibrators with accurate 
independent distances. The slope is given by the regression with errors in linewidths to 
the 151 cluster galaxies. The separate components of the plot are more easily seen in the 
figures that follow. 



the calibration slope does not depend on the galaxy luminosity function. If 
there was evidence of a slope change, a slightly more complicated analysis 
involving a non-linear relationship would have been necessary. 

3.3. ABSOLUTE CALIBRATION AT I BAND 

Curently there are 17 galaxies with distances determined through observa- 
tions of cepheid variable stars, mostly from observations with the Hubble 
Space Telescope (Freedman et al. 1997, Sandage et al. 1996, Tanvir et al. 
1995). Two of the systems, NGC 2366 and NGC 3109, are fainter than the 
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Figure 3. Luminosity-linewidth relation for distance calibrators. 
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Figure 4- Luminosity-linewidth relation for Ursa Major Cluster. 



template cutoff so will be ignored. One more galaxy is added as a calibra- 
tor, NGC 4258, which has a distance from the geometry inferred for the 
circum-nuclear masers (Miyoshi et al. 1995). Hence 16 calibrators are used. 

It would be improper to do a regression on the calibrator relationship 
because in no way do they provide a complete sample. It can only be as- 
sumed that the calibrators are drawn from a similar distribution as the 
template objects, perhaps with magnitude as a selection criterion but not 
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Figure 5. Luminosity-linewidth relation for Fornax Cluster. 
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Figure 6. Luminosity-linewidth relation for Pisces filament. 



linewidth. Effectively, each of the 16 calibrators provides a separate zero- 
point offset. The least-squares average provides the optimum fit. The final 
result is shown in Figure 2 where the / band luminosity-linewidth relation 
is shown for the 16 calibrators and the 151 cluster template galaxies shifted 
to the absolute magnitude scale of the calibrators (the 2 rejected galaxies 
are also plotted). Figures 3-8 present the same material but separated to 
distinguish the fits to the calibrators and the individual clusters. 
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Figure 7. Luminosity-linewidth relation for Abell 1367 Cluster. 
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Figure 8. Luminosity-linewidth relation for Coma Cluster. 
3.4. THE B, R AND K' RELATIONS 

Less complete information is available at other bands than /. However, 
inter-comparisons are valuable because of the potential problem with ob- 
scuration. Information is available at B and R for the calibrators, all the 
galaxies in the Ursa Major sample, most of those in Coma, and for most 
in the part of the Pisces region at 00 h 49 m < a < 0l h 32 m (Pierce & Tully 
1998). Material is available at K' for the same Ursa Major and Pisces 
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Figure 9. Luminosity-linewidth relations in B, R, I, K' after corrections for inclination. 
Filled circles: Ursa Major; open circles: Pisces. Straight lines are regressions with errors 
in linewidths. 
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galaxies (Tully et al. 1996, 1998). The luminosity-linewidth relations are 
illustrated in Figure 9. 

The magnitude scatter is essentially the same at R, I, K' and ~ 20% 
worse at B. Obscuration corrections diminish toward the infrared until 
they are tiny at K'. However, sky background contamination causes degra- 
dation toward the infrared from the favorable situation at R to the poor 
situation at K' where one looses almost 2 scalelengths to the sky compared 
with an R exposure of the same duration. The correlations are seen to 
steepen toward the infrared. However, this steepening is less extreme than 
had been seen in the past because of the strong luminosity dependence 
of the reddening corrections that are now applied. The biggest corrections 
are made to the most luminous galaxies in the bluest bands. Hence the 
corrected relations at shorter wavelengths are steepened toward the slopes 
of the almost-reddening- free infrared relations. As shown in Tully et al. 
(1998), only a weak color dependency on luminosity remains after redden- 
ing is taken into account. Slopes at B, R, I, K' are -7.8, -8.0, -8.2, and -8.7, 
respectively, with the correlation against the same linewidth information. 
These slopes are based on the linewidth regression which is appropriate 
for bias-free distance determinations but not the slopes that one wants to 
give a physical interpretation. The true nature of the correlation is char- 
acterized better by a maximum likelihood fit. An approximation to that 
is a double regression, which at / gives the slope -7.9. These fits indicate 
infrared convergence toward L oc W n where n = 3.4 ± 0.1. 

4. Summary of Results 

The template relation with the zero-point given by 16 galaxies can be used 
to determine distances to any other galaxy, with the proviso that the procee- 
dure must fail if the target is intrinsically less luminous than Mj ,l,k = — 18 m , 
the faintness limit of the calibration. If the problem is to measure Ho, this 
faint limit is of no concern because targets of interest are beyond the Local 
Supercluster where peculiar velocities are expected to be a small fraction 
of expansion velocities. 

It would be possible to apply the calibration to measure distances to 
hundreds of galaxies. For the moment, with the interest of maintaining as 
homogeneous a set of measurements as possible, the Ho determination will 
be based on the 5 clusters that went into the template plus 7 other clusters 
each with of order a dozen measures. The results are presented in Table 1 
and Figure 10. The table provides (col. 2) the number of measures in the 
cluster, (col. 3) the rms scatter about the template relation, (col. 4/5) the 
distance modulus/distance of the cluster, (col. 6) the velocity of the cluster 
in the CMB frame as given by Giovanelli et al. (19976), and (col. 7) the 
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measure of Ho from the cluster. The velocity given to the Pisces filament 
is the average of the values for the three main sub-condensations. 



TABLE 1. Five Template Clusters and Seven More 



Cluster 


No. 


RMS 


Modulus 


Distance 


V^m6 


Ho 






(mag) 


(mag) 


(Mpc) 


(km/s) 


(km/s/Mpc) 


Fornax 


16 


0.50 


31.21 


17.5 


1321 


76 


Ursa Major 


38 


0.41 


31.34 


18.5 


1101 


59 


Pisces Filament 


46 


0.31 


33.90 


60.3 


(4779) 


79 


Abell 1367 


23 


0.41 


34.64 


84.8 


6735 


79 


Coma 


28 


0.33 


34.67 


85.9 


7185 


84 


Antlia 


11 


0.27 


32.78 


35.9 


3120 


87 


Centaurus 30 


13 


0.52 


32.94 


38.9 


3322 


86 


Pegasus 


12 


0.37 


33.36 


46.9 


3519 


75 


Hydra I 


11 


0.35 


33.84 


58.6 


4075 


70 


Cancer 


16 


0.34 


33.89 


60.0 


4939 


83 


Abell 400 


9 


0.24 


34.91 


96.1 


6934 


72 


Abell 2634 


16 


0.32 


35.23 


111.0 


7776 


70 


Weighted average 












77 



The error bars in Fig. 10 contain distance and velocity components. The 
errors associated with distance depend directly on the rms dispersion in a 
cluster and inversely with the square root of the number of galaxies in the 
cluster sample. The error associated with velocity streaming is taken to be 
300 km s _1 . The velocity component to the error is totally dominant inside 
2000 km s _1 . The statistical errors in distance become the dominant factor 
beyond ~ 6000 km s _1 . The symbols in Fig. 10 differ for different regions 
of the sky. There is a hint of systematics: for example the filled circles lie 
above the open circles. For the present purposes, the best estimate of Ho is 
derived by taking an average of log Ho values with weights proportional to 
the inverse square of the error bars that are plotted. The result is Ho = 77 ± 
4 km s _1 Mpc -1 . The error is the 95% probability statistical uncertainty. 
It is small because it is based on a template of 151 galaxies, a zero-point 
fixed by 16 galaxies, and application to 12 clusters distributed around the 
sky and out to 8,000 km s" 1 . 

This result is somewhat higher than 'interm' value of Ho ~ 73 reported 
by the HST Key Project team (Mould et al. 1997) or the value of H = 
69 ± 5 found by Giovanelli et al. 1997a based on similar applications of 
luminosity-linewidth correlations. The latter value is outside the bound of 
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Figure 10. Individual estimates of Ho as a function of systemic velocity. Errors are 
a convolution of the statistical errors in distance and an uncertainty of 300 km s _1 in 
velocities. Symbols vary with location on the sky: filled circles: north celestial and north 
galactic; open circles: north celestial and south galactic; filled squares: south celestial and 
north galactic; open squares: south celestial and south galactic. 



the statistical error found in this paper. It must be attributed to a small 
systematic difference that remains to be identified. 

The value of Ho = 77 found here is about 13% lower than the value 
found by the same author and method in the past. There has been a 17% 
decrease due to the revision of the luminosity-linewidth zero-point as a con- 
sequence of the increase from 4 to 16 in the number of calibrator galaxies 
with distances determined through the cepheid period-luminosity relation. 
There has been an 8% increase as a consequence of the revised reddening 
corrections now being applied. A 5% decrease has come about with the 
introduction of material on more clusters around the sky at distances well 
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beyond the Local Supercluster. These changes are a sobering illustration 
of random and systematic errors. The shift associated with the improved 
cepheid calibration is comparable to the rms dispersion of the luminosity- 
linewidth relations at R, I, K' bands. Either by statistical fluke or an un- 
known systematic, the original 4 calibrators are drawn from the faint side 
of the intrinsic correlation. 

In conclusion, there has been a decrease of ~ one standard devia- 
tion from the value of Ho measured previously by this author with the 
luminosity-linewidth method, to Ho = 77 ±4 km s" 1 Mpc -1 . The greatest 
perceived problems of old have been addressed: there are now many more 
zero-point calibrators, the template is much more extensive and complete, 
reddening corrections are under better control, and the method is applied 
to more targets distributed around the sky. Formal statistical errors are cut 
in half. Uncertainties are now dominated by potential, unidentified system- 
atics. 
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